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Abstract. For a nonsingular integer matrix A, we study the growth of the 
order of A modulo N. We say that a matrix is exceptional if it is diagonalizable, 
and a power of the matrix has all eigenvalues equal to powers of a single rational 
integer, or all eigenvalues are powers of a single unit in a real quadratic field. 

For exceptional matrices, it is easily seen that there are arbitrarily large 
values of A'^ for which the order of A modulo A'^ is logarithmically small. In 
contrast, we show that if the matrix is not exceptional, then the order of A 
modulo N goes to infinity faster than any constant multiple of log N. 



1. Introduction 

Let ^ be a c? X d nonsingular integer matrix, and iV > 1 an integer. The order, 
or period, of A modulo N is defined as the least integer A: > 1 such that A'^ = I 
mod N, where I denotes the identity matrix. If A is not invertible modulo N then 
we set ord(^, N) — oo . In this note we study the minimal growth of ord(^, N) as 
N ^ oo. 

If A is of finite order (globally), that is A^ ~ I for some r > 1, then clearly 
oid{A,N) < r is bounded. If A is of infinite order, then oi-d{A,N) — > oo as 
N oo. Moreover, in this case it is easy to see that ord(v4, N) grows at least 
logarithmically with N, in fact if no eigenvalue of A is a root of unity then: 

OTd{A,N) > — logiV + 0(l) 
VA 

where tja '■— X]|a |>i ^"^S I'^il' sum over all eigenvalues {A^} of A which lie 
outside the unit circle {ija is the entropy of the endomorphism of the torus K'^/Z'^ 
induced by A, or the logarithmic Mahler measure of the characteristic polynomial 
of A, and the condition that no eigenvalue of A is a root of unity is equivalent to 
ergodicity of the toral endomorphism). 

There are cases when the growth of ord(A, N) is indeed no faster than logarith- 
mic. For instance if we take d = 1, and A = (a) where a > 1 is an integer, and 



logiVfc 
log a 



ord(A, Nk) = k 
and so 

n^ r • .ord(A,iV) 1 

(1) limmf ^ = < oo 

logiV loga 

in this case. 
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The same behaviour occurs in the case of 2 x 2 unimodular matrices A G SL2 (Z) 
which are hyperbohc, that is A has a pair of distinct real eigenvalues A > 1 > A~^. 
Then 

(2) hminf-^fY^^A-- 

logiV log A ?u 

See e.g. IETT2I ^. 

These cases turn out to be subsumed by the following definition: We say that 
A is exceptional if it is of finite order or if it is diagonalizable and a power A^ of A 
satisfies one of the following: 

(1) The eigenvalues of A^ are all a power of a single rational integer a > 1; 

(2) The eigenvalues of A^ are all a power of a single unit A 7^ ±1 of a real 
guadratic field. 

We will see that if A is exceptional, then there is some c > and arbitrarily 
large integers N for which ord{A, N) < clog TV. 
Our main finding in this note is 

Theorem 1. If A ^ Mat(i(Z) is not exceptional then 

ord(A, N) 

> 00 

logA^ 

as N ^ 00. 

A special case is that of diagonal matrices, e.g. ^ = ^ . In that case The- 

orem^says that ord(a, 6; A'^)/ logiV —> 00 if a, 6 are multiplicatively independent, 
in contrast with 

Theorem ^ is in fact equivalent to a subexponential bound on the greatest com- 
mon divisor gcd(A" — /) of the matrix entries of A" — I. We shall derive it from 

Theorem 2. If A E Matci(Z) is not exceptional then for all e > 

gcd(A" - /) < exp(en) 

if n is sufficiently large. 

In the special case of a diagonal matrix such a,s A — ^ , we have gcd(A" — 

/) = gcd(a" — 1,6" — 1). In |BCZ | it is shown that if a,b are multiplicatively 
independent then for all e > 0, 

(3) gcd(a" - 1, 6" - 1) < exp(en) 

for n sufficiently large, giving Theorem [51 in that case. To prove Theorem 12 in 
general, we will use a version of (|3J for 5-units in a general number field |CZ| . 

We note that Theorem |21 establishes upper bounds on gcd(A" — /). As for lower 
bounds, it is conjectured in |AR) that if A has a pair of multiplicatively independent 
eigenvalues then liminf gcd(A" — /) < 00. 

Motivation: A natural object of study for number theorists, the periods of 
toral automorphisms were also investigated by a number of physicists and mathe- 
maticians interested in classical and quantum dynamics, see e.g. |HBI IR1 FPFI . One 



special case of this appeared as a problem in the 54-th W.L. Putnam Mathematical Com- 
petition, 1994, see Eil Pages 82, 242]). 
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reason for our own interest also lies in the quantum dynamics of toral automor- 
phisms: It has recently been shown that any ergodic automorphism A G SL2(Z) 
of the 2-torus admits "quantum limits" different from Lebesgue measure [FNB" , if 
one does not take into account the hidden symmetries ("Hecke operators") found 
in |KR1| . The key behind the constructions of these measures is the existence of 
values of N satisfying (j^J, that is ord{A, N) ^ 2\ogN /rjA- A higher-dimensional 
version of this would involve taking ergodic symplectic automorphisms A G Sp2g (Z) 
of the 2f;-dimensional torus. Theorem^gives one obstruction to extending the con- 
struction of jFNBj to the higher-dimensional case. 

2. Proof of Theorem [21 

Assume that for a certain positive e and all integers n in a certain infinite se- 
quence TV C N we have 

(4) gcd(A" -/) > exp(en). 

We shall prove that A is "exceptional" , in the sense of the above definition. 

We let fc C Q be the splitting field for the characteristic polynomial of A, so we 
may put A in Jordan form over fc, namely, we may write 

A = PBP-\ 

where P is an invertible d x d matrix over fc and B is in Jordan canonical form. 

For later reference we introduce a little notation related to the field fc. 

We let M (resp. Mq) denote the set of (resp. finite) places of fc. We shall 
normalize all the absolute values with respect to fc, i.e. in such a way that the 
product formula Il^gj^j = 1 holds for x G fc*, and the absolute logarithmic 
Weil height reads h{x) = logmaxjl, jxl^}. We also let S" be a finite set of 
places of fc including the archimedean ones and we denote by Og the group of 
S-units in fc*, namely those elements x € k such that |a:|^ = 1 for all ^ ^ S. 

Note that — / = P~^{A^ — I)P; since the entries of P and its inverse are 
fixed independently of n, hence have bounded denominators as n varies, this formula 
shows that the entries of — / have a "big" g.c.d., in the sense of ideals of fc, for 
n G Af. Since the entries of — I are algebraic integers, not necessarily rational, 
to express their g.c.d. we shall use the formula-definition 

loggcd(B" - /) := V log- max - 

where log~(a;) := — min(0, loga;); this is a nonincreasing nonnegative function of 
x>0. 

Note that this definition agrees with the usual notion in case B has rational 
integer entries. From Q and the above formula _B" — / = P^^(yl" — I)P we 
immediately deduce that 

(5) log- max - /),y|„ > for large n G AA . 

— ' ij 2 

Me Mo 

In fact, each entry of i?" — / is a linear combination of entries of A" — I with 
coefficients having bounded denominators, whence |(-B" — < c^inax^s \ {A" — 

I)rs\fj., where are positive numbers independent of n such that = 1 for all but 
finitely many fj, G A/. This proves 
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We start by showing that B must be necessarily diagonaL In fact, if not some 
block of B would contain on the diagonal a 2 x 2 matrix of the form 

'\ I' 
.0 A, 

where A is an (algebraic integer) eigenvalue of A. Hence _B" — / would contain 
among its entries the numbers A" — 1 and A"^^n. Then, for every /x G Mq, we 
would have 

max|(B"-/),,|^ >max(|A"-l|,„|A"-in|,J > 

whence log" max^j \ {B" — I)ij\f_i < log^ = — log \n\^. In conclusion, 

V log" max|(S" - < V -log|n|^ = logn 
— — 

the last equality holding because of the product formula. However this contradicts 
© for all large n e A/" and this contradiction proves that B is diagonal. 

Therefore from now on we assume that _B is a diagonal matrix formed with the 
eigenvalues Ai, . . . , A^ of A, each counted with the suitable multiplicity. 

Another case now occurs when there exist two multiplicatively independent 
eigenvalues, denoted a,/3. Now, from (0) we get, for large n E JV, 

(6) log" max(|a" - 1^, - 1^) > ^ log" max - > . 

We are then in position to apply (after a little change of notation) the following 
fact from |UZj . stated as Proposition 2 therein: 

Proposition 3 (Proposition 2 of jCZ| ). Let 6 > 0. All but finitely many solutions 



(w, v) G (Og)'^ to the inequality 



log max{|u — |w — l|p} > (5 • max{ft,(u), /i(w)} 

fie Mo 

satisfy one of finitely many relations u"'v^ — 1, where a,b G Z are not both zero. 

Actually, Prop. 2 in [CZ] is a little stronger, since the summation is over all 
H G M rather than the finite fj, € Mq and since it also asserts that the relevant 
pairs (a, b) may be computed in terms of 6. 

We apply this fact with u — a'"', v — (3"^ and S containing the finite set of places 
of k which are nontrivial on a or /3; note that © implies the inequality of the 
proposition, with 6 = e / {2 max{h{a) , h{f3))) . We conclude that, for an infinity of 
n e J\f, a same nontrivial relation a""/?''" — 1 holds, contradicting the multiplicative 
independence of a, /?. 

Therefore we are left with the case when all pairs of eigenvalues are multiplica- 
tively dependent. This means that they generate in k* a subgroup F of rank < 1. 

If the rank is zero all the eigenvalues Ai are roots of unity, so the matrix A has 
finite order and thus it is exceptional. Hence let us assume from now on that the 
rank is 1. Let then A S F be a generator of the free part of F (it exists by basic 
theory) . Then, for suitable roots of unity (^i , . . . , Cd and rational integers ai , . . . , ad 
we may write 

(7) A, = C»A"', i = l,...,d. 
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Necessarily the Ci lie in k. 

Let a be an automorphism of k. Then a fixes the set of eigenvalues, since A is a 
matrix defined over Q; hence a fixes the above group F. Let r be the order of the 
torsion in T, so the subgroup [r]T of r-th powers in F is cyclic, generated by A*". 
(Note that automatically = 1 in Q). Then cr must send A"" to another generator 
of [r]F, whence 

a{\Y = A^'-. 

Therefore in particular A*" is at most quadratic over Q (in fact, recall that k/Q is 
normal) . 

Let us first assume that A*" is rational. Raising the equations Q to the power 
2r, we see that the eigenvalues A^'' of the matrix A^^ are positive rationals; since 
they are algebraic integers, they are therefore positive rational integers. Since they 
are pairwise multiplicatively dependent they are powers of a same positive integer 
(which can be taken A^^''). We thus fall in another of the exceptional situations. 

The last case occurs when A'' is a quadratic irrational. Then some automorphism 
a must send it to its inverse A^*". As before, we may raise equations Q to the r-th 
power to find A[ = A''"'. Therefore (t(A[) = A"*^. Since the Xi are algebraic integers, 
the same is true for the A^'^, and hence we find that all the eigenvalues of are 
units (some of them possibly equal to ±1) in a same quadratic field. 

This concludes the proof. 



3. Proof of Theorem [I] 
The following Lemma shows that Theorems ^ and |21 arc in fact equivalent: 

Lemma 4. Let A be a nonsingular integer matrix of infinite order. Then the 
following are equivalent: 

(1) For all e > 0, we have gcd(yl" — /) < exp(en) if n is sufficiently large; 

(2) ord(A,iV)/logiV -> oo. 

Proof. Assume that gcd(^" — I) < exp(en) for all e > 0. Fix e > 0. Take n = 
ord(v4, N) and note that N divides all the matrix entries of ^°''d(A,Af) _ j Since 
A does not have finite order and thus ord(yl; N) ^ oo as — > oo, we have for A^ 
sufficiently large that 

N < gcd(yl°'''^(^'^) - /) < exp(eord(yl, A^)) 

Thus 

log A' < eord(^. A') . 

Since this holds for all e > we find ord{A, N) / log N co. 

Conversely, suppose that there is some p > and an infinite sequence of integers 
Af so that gcd(A" — /) > exp(pn) for all n G JV. Then for the sequence A'n :— 
gcd(yl" — I), n ^ JV (which is infinite since A^„ > exp(pn)) we have 

ord(^, Nn)<n< loggcd(yl" - I)/p = log N^/p 

and thus lim inf ord(A, N)/ logN < oo. □ 
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4. Comments 

It is readily seen that exceptional cases do in fact occur, and that they give rise 
to powers A" such that gcd(v4" — /) is exponentially large, and hence to arbitrarily 
large integers N for which ord(^, N) is logarithmically small. The last case of the 
eigenvalues in a quadratic field of course requires that the irrational ones occur in 
conjugate pairs, since A is defined over Q, and that the determinant of A is ±1. 
Examples of such integer matrices can be produced from the action of a fixed such 
2x2 hyperbolic matrix Aq g S'L2(Z) on tensor powers, or from Aq iSi cr where a is 
a permutation matrix. 

To see that the exceptional cases lead to exponentially large gcd, consider first 
the case that a power of A has all eigenvalues a power of a single integer a > 1. As 
we have seen in the course of proof of Theorem[21 replacing a matrix by a conjugate 
(over Q) does not change the asymptotic behaviour. Thus we may assume that is 
diagonal with eigenvalues a™"^, . . . , a™''. Then clearly ord{A^, N) < ord(a, N) and 
taking Nn := a" — 1 gives ord(a, 7V„) = n ^ log Nn/ a. Thus we find ord{A, iV„) < 
r\ogNn/a. 

Now assume that a power A''' of A has all its eigenvalues a power of a single 
unit A > 1 in a real quadratic field K. Then for some matrix P with entries in K , 
we have A^ = PBP'^ with B diagonal with eigenvalues A°^, . . . , A"'', where are 
integers which sum to zero. 

Since P is only determined up to a scalar multiple, we may, after multiplying P 
by an algebraic integer of K , assume that P has entries in the ring of integers Ok 
of K, and then P~^ = gj^j^-P"'* where P'"' also has entries in Ok- 

The entries of A^'' — I are thus Ox-linear combinations of (A"''^ ^1)/ det(P). 
We now note that 

A-''' - 1 = -A-'^(A'= - 1) 

and thus the entries of A^^ — I are all Ox-linear combinations of (A'"'''^ — 1)/ det(P), 
which are in turn Ox-multiples of (A'^ ^1)/ dct(F). In particular, gcd(A'"'^ — /), 
which is a Z-linear combination of the entries of A^'^ — I, can be written as 

with "fk e Ok- 

Now taking norms from to Q we see 



Since 7^. ^ 0, we have |Nx/Q(7fe)| > 1 and thus 

which gives | gcd(A''''' — /)| ^ A*^/^, namely exponential growth. 
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